Antiferromagnetic Ising model in small- world networks 
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The antiferromagnetic Ising model in small-world networks generated from two-dimensional regu- 
lar lattices has been studied. The disorder introduced by long-range connections causes frustration, 
which gives rise to a spin-glass phase at low temperature. Monte Carlo simulations have been car- 
ried out to study the paramagnetic to spin-glass transition, as a function of the rewiring probability 
p, which measures the disorder strength. The transition temperature T c goes down for increasing 
disorder, and saturates to a value T c w 1.7 J for p > 0.4, J being the antiferromagnetic coupling. 
For small p and at low temperature, the energy increases linearly with p. In the strong-disorder 
limit p — * 1, this model is equivalent to a short-range ±J spin glass in random networks. 

PACS numbers: 64.60.De, 05.50.+q, 75.10.Nr, 89.75.Hc 



I. INTRODUCTION 

In the last few years, there has been a surge of interest 
in modeling complex systems as networks or graphs, with 
nodes representing typical system units and edges play- 
ing the role of interactions between connected pairs of 
units. Thus, complex networks have been used to model 
several types of real-life systems (social, economic, bi- 
ological, technological), and to study various processes 
taking place on them [T|, Q, HI, EL H| • In this context, some 
models of networks have been designed to explain empir- 
ical data in several fields, as is the case of the so-called 
small- world networks, introduced by Watts and Strogatz 
in 1998 @. 

These small-world networks are well suited to study 
systems with underlying topological structure ranging 
from regular lattices to random graphs @, H[ , by chang- 
ing a single parameter Q. They are based on a regular 
lattice, in which a fraction p of the links between nearest- 
neighbor sites are replaced by new random connections, 
creating long-range "shortcuts" [E E|- In the networks 
so generated one has at the same time a local neigh- 
borhood (as in regular lattices) and some global prop- 
erties of random graphs, such as a small average topo- 
logical distance between pairs of nodes. These networks 
are suitable to study different kinds of physical systems, 
as neural networks [Tol ] and man-made communication 
and transportation systems [E [Til [l2j • The importance 
of a short global length scale has been emphasized for 
several statistical physical problems on small-world net- 
works. Amongthese problems, one finds the spread of 
infections [IE [3] , signal propagation lrllT^.[T6l , random 
spreading of information U/A Il8l. [l9l 120 . l2l| , as well as 
site and bond percolation [14 , |22| . [23[ . 

Cooperative phenomena in this kind of networks are 
expected to display unusual characteristics, associated 
to their peculiar topology [13, HE HE Thus, a 

paramagnetic-ferromagnetic phase transition of mean- 
field type was found for the Ising model on small- 
world networks derived from one-dimensional (ID) lat- 
tices [H], HE HIl- This phase transition occurs for any 



value of the rewiring probability p > 0, and the transition 
temperature T c increases as p is raised. A similar mean- 
field-type phase transition was found in small-world net- 
works generated from 2D and 3D regular lattices [IE HT[ , 
as well as for the XY model in networks generated from 
one-dimensional chains 32]. In recent years, the Ising 
model has been thoroughly studied in complex networks, 
such as the so-called scale-free networks, where several 
unusual features were observed [H, HE HE HE] • 

Here we study the antiferromagnetic (AFM) Ising 
model in small-world networks generated by rewiring a 
2D square lattice. One expects that the AFM ordering 
present in the regular lattice at low temperature will be 
lost when random connections are introduced, for an in- 
creasing number of bonds will be frustrated as p rises. 
In particular, this model includes the two basic ingredi- 
ents necessary to have a spin glass (SG), namely, disorder 
and frustration. The former appears due to the random 
long-range connections introduced in the rewiring pro- 
cess, and the latter because half of these rewired links 
connect sites located in the same sublattice of the start- 
ing regular lattice. 

In some spin-glass models, such as the Sherrington- 
Kirkpatrick model, all spins are mutually connected 
[37L HE]. An intermediate step between these globally 
connected networks and finite-dimensional models con- 
sists in studying sp in g lasses on random graphs with finite 
(low) connectivity [39j, [40j, |4lJ, |42j . A further step between 
random graphs with finite mean connectivity and regu- 
lar lattices is provided by small- world networks, where 
one can modify the degree of disorder by changing the 
rewiring probability p. Then, for the AFM Ising model 
on small-world networks, we expect to find features close 
to those of short-range spin-glass systems. In this line, a 
spin-glass phase has been recently found and character- 
ized for the AFM Ising model in scale- free networks [43[ . 
In this paper, we employ Monte Carlo (MC) simulations 
to study the paramagnetic to spin-glass phase transition 
occurring in small-world networks. Apart from temper- 
ature and system size, another variable is the rewiring 
probability, which controls the degree of disorder, and 
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allows us to interpolate from a paramagnetic- AFM tran- 
sition at p = to a paramagnetic-SG transition in a 
random graph at p = 1 . 

The paper is organized as follows. In Sec. II we de- 
scribe the networks and the computational method em- 
ployed here. In Sec. Ill we give results for the heat ca- 
pacity, energy, and spin correlation, as derived from MC 
simulations. In Sec. IV we present and discuss the overlap 
parameter, transition temperature, and absence of long- 
range ordering. The paper closes with the conclusions in 
Sec.V. 



II. MODEL AND METHOD 



We consider the Hamiltonian: 

H = JijSiSj , 



(1) 



where Si — ±1 (i = l,...,N), and the coupling matrix 
Jij is given by 

7 _ J J(> 0) if i and j are connected, 
J y = i otherwise. [ ' 



This means that each edge in the network is an AFM 
interaction between spins on the two linked nodes. Note 
that, contrary to the usually studied models for spin 
glasses, in this model all couplings are antiferromagnetic. 
This model with AFM couplings can be mapped onto one 
in which all unrewired bonds are ferromagnetic (FM), 
and the rewired links are 50% FM and 50% AFM. In the 
limit p — > 0, this mapping is the well-known correspon- 
dence between AFM and FM Ising models on bipartite 
lattices (44[. In the limit p — > 1, our AFM Ising model 
is equivalent to a spin-glass model on a random graph of 
mean connectivity (k) = 4, with 50% AFM and 50% FM 
bonds. 

Small-world networks have been built up according to 
the model of Watts and Strogatz [f| i.e., we consid- 
ered in turn each of the bonds in the starting 2D lattice 
and replaced it with a given probability p by a new con- 
nection. In this rewiring process, one end of the selected 
bond is changed to a new node chosen at random in the 
whole network. We impose the conditions: (i) no two 
nodes can have more than one bond connecting them, 
and (ii) no node can be connected by a link to itself. 
With this procedure we obtained networks where more 
than 99.9% of the sites were connected in a single com- 
ponent. Moreover, this rewiring method keeps constant 
the total number of links in the rewired networks, and 
the average connectivity (k) coincides with z = 4. This 
allows us to study the effect of disorder upon the physi- 
cal properties of the model, without changing the mean 
connectivity. 

We note that other ways of generating small- world net- 
works from regular lattices have been proposed [13, [45| . 



In particular, instead of rewiring each bond with prob- 
ability p, one can add shortcuts between pairs of sites 
taken at random, without removing bonds from the reg- 
ular lattice. This procedure turns out to be more conve- 
nient for analytical calculations, but does not keep con- 
stant the mean connectivity (fc), which in this case in- 
creases with p. Spin glasses on such small-world net- 
works, generated from a one-dimensional ring, have been 
studied earlier by replica symmetry breaking [46| and 
transfer matrix analysis [4?| . 

From the 2D square lattice, we generated small- world 
networks of different sizes. The largest networks em- 
ployed here included 200 x 200 nodes. Periodic boundary 
conditions were assumed. For a given network, we carried 
out Monte Carlo simulations at several temperatures, 
sampling the spin configuration space by the Metropolis 
update algorithm (48|, and using a simulated annealing 
procedure. Several variables characterizing the consid- 
ered model have been calculated and averaged for differ- 
ent values of p, T, and system size N. In general, we have 
considered 300 networks for each rewiring probability p, 
but we used 1000 networks to determine accurately the 
transition temperature from paramagnetic to SG phase. 
In the following, we will use the notation (...) to indicate 
a thermal average for a network, and [...] for an average 
over networks with a given degree of disorder p. 



III. THERMODYNAMIC OBSERVABLES 

The heat capacity per site, c v , was obtained from the 
energy fluctuations AE at a given temperature, by using 
the expression 



[(A#) 2 
NT 2 



(3) 



where (AE) 2 = (E 2 ) - (E) 2 . We have checked that the 
results coincide within numerical noise with those derived 
by calculating c v as [d(E] /dT)/N . Note that we take the 
Boltzmann constant kB = 1- 

The temperature dependence of c v is displayed in Fig. [I] 
for saveral values of the rewiring probability p and for 
networks built up from a 80 x 80 2D lattice. For increas- 
ing p, one observes two main features: the maximum of 
c v shifts to lower T and the peak broadens appreciably. 
This broadening agrees with the behavior expected for 
systems with increasing disorder, similarly to that found 
for the FM Ising model in these networks [30]. However, 
in the AFM model, the shift of the peak to lower tem- 
perature suggests a phase transition with a temperature 
T c that decreases as p is raised, contrary to the FM case, 
where an increase in T c with p was observed. This differ- 
ence between both Ising models on small- world networks 
occurs in addition to the nature of the transition itself, 
which in the FM case is a paramagnetic-ferromagnetic 
transition vs a paramagnetic-SG transition in the AFM 
model (see below). A decrease in T c for the AFM Ising 
model in this kind of networks was also suggested in Ref. 
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FIG. 1: Heat capacity per site c„ vs temperature for small- 
world networks generated from a 2D lattice of size 80x80. The 
plotted curves correspond to different values of the rewiring 
probability p. From top to bottom: p — 0.1, 0.3, 0.5, and 1. 



149 from the behavior of the heat capacity for several val- 
ues of p. 

The increase in disorder as p is raised is accompanied 
by an increase in frustration of the links at low temper- 
atures. This can be quantified by the low-temperature 
energy of the system, which will rise as the rewiring 
probability is raised. To obtain insight into this energy 
change for p near zero (low disorder), let us remember 
that the square lattice is bipartite, in the sense that one 
can define two alternating sublattices, say A and B, so 
that neighbors of each node in sublattice A belong to 
sublattice B, and vice versa. In the rewiring process, 
one introduces links between nodes in the same sublat- 
tice, and the resulting networks are no longer bipartite. 
However, for small rewiring probability p, we can still 
speak about two sublattices, with some "wrong" connec- 
tions. Since each link is rewired with probability p, each 
connection in the starting regular lattice will be trans- 
formed into a wrong connection (of types A-A or B-B) 
with probability p/2. The remaining links are of A-B 
type, and the number of wrong connections is on aver- 
age zNp/A. Then, for small p, the lowest energy can 
be approximated by E m = —zNJ(l — p)/2, under the as- 
sumption that the AFM long-range ordering of the square 
lattice is still preserved. For z = 4, we have an energy 
per node: e m = —2(1 — p) J. Note that for finite p the 
low-temperature long-range ordering in fact decays due 
to the appearance of domains driven by the rewired con- 
nections (see below), but the AFM ordering is a good 
reference to obtain insight into the energy change as a 
function of rewiring probability p. 

We now turn to the results for the minimum energy 
reached in our simulations for different p values, which 
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FIG. 2: Minimum energy per site obtained in our simulations 
for the AFM Ising model on small-world networks with N = 
10 4 nodes. The dashed line corresponds to e m = (—2 + 2p) J, 
as explained in the text. The dotted line indicates the ground- 
state energy obtained in Ref. |4l| for a spin glass on random 
networks with (fe) = 4. The solid line is a guide to the eye. 



are shown in Fig. [2] For rising p, E^ increases from the 
value corresponding to AFM ordering in the regular lat- 
tice, e m = — z J/2. The dashed line in Fig. [2] displays 
the behavior expected for small p, in the case of a strict 
AFM ordering on the underlaying lattice. This estima- 
tion is close to the minimum energy obtained in our sim- 
ulations for p < 0.15. For larger p values, it departs 
appreciably from the results of the simulations, and e m 
lies below the dashed line. In the limit p — 1 we find a 
value e m = — 1.444(2)J. In this limit, our small-world 
networks are very close to random networks with a Pois- 
son distribution of connectivities, but are not identical to 
the latter because of the restriction that no nodes have 
zero links, imposed in the rewiring process flU [24| ■ For 
a ± J Ising spin glass on random graphs with a Poisson 



distribution of connectivities and (k) = 4, Boettcher 41] 
found a ground-state energy e m = —1.431(1) J by using 
extremal optimization. This value is plotted in Fig. [2] as 
a dotted line close to p — 1 , and is near the minimum en- 
ergy we found for the small- world networks in this limit. 
For a more direct comparison with random networks, we 
have carried out some simulations for small-world net- 
works with p = 1, where we allowed the presence of 
isolated sites (with connectivity k = 0). For the AFM 
Ising model in these networks, we found a minimum en- 
ergy e m = — 1.438(2)J, between those of our standard 
networks (with minimum connectivity k — 1) and ran- 
dom networks in Ref. l4ll Note that our error bar in e m 
corresponds to a standard deviation in the distribution 
of minimum energy obtained for different networks. We 
emphasize that the energy e m found here for each value 
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FIG. 3: Absolute value of the spin correlation function vs 
distance in small-world networks with N = 10 4 nodes, at 
temperature T = 1.5J. The dimensionless distance between 
nodes, r — d/do, is measured on the starting regular lattice. 
Data are shown for several values of p, as indicated by the 
labels. 



of p is an upper limit for the lowest energy of the system. 

Even though the random connections present in small- 
world networks introduce disorder in the starting regular 
lattice, these networks still keep memory of the original 
bipartite lattice, but the actual meaning of the partition 
in two sublattices is gradually reduced as p rises. This 
can be measured by the number of wrong links, which 
amounts to a fraction p/2 of the total number of links, as 
indicated above. In the limit p = 1, half of the links con- 
nect sites in the same original sublattice, and the mem- 
ory of the partition in sublattices has completely disap- 
peared. One can visualize the loss of AFM ordering on 
the 2D lattice, by plotting the spin correlation vs distance 
for several values of p. We define £ as 



CM - [(SiSj) 



(4) 



where the subscript r indicates that the average is taken 
for the ensemble of pairs of sites at distance r. Note 
that r = d/do refers here to the dimensionless distance 
between sites in the starting regular lattice, not to the 
actual topological distance or minimum number of links 
between nodes in the rewired networks (do is the dis- 
tance between nearest neighbors). The correlation £(r) 
is shown in Fig. [3] for several values of the rewiring prob- 
ability p, at temperature T = 1.5 J. This temperature 
is below the critical temperature T c of the paramagnetic- 
SG transition for all values of p (see below). As expected, 
£(r) decreases faster for larger p, and vanishes at p = 1 
for any distance r > 1. In general, after a short tran- 
sient for small r, we find an exponential decrease of the 
spin correlation with the distance. This indicates that, 



in spite of the disorder present in the networks for p > 0, 
there remains some degree of short-range AFM ordering 
on the starting regular lattice, which is totally lost in the 
limit p — ► 1. 



IV. SPIN-GLASS BEHAVIOR 

A. Overlap parameter 

As is usual in the study of spin glasses, we now con- 
sider two copies of the same network, with a given real- 
ization of the disorder, and study the evolution of both 
spin systems with different initial values of the spins and 
different random numbers for generating the spin flips 
[Hoi . [5lT | . It is particularly relevant the overlap q between 
the two copies, defined as 



(1) 5(2) 



(5) 



Ob- 



where the superscripts (1) and (2) denote the copies, 
viously, q is defined in the interval [—1,1]. 

We have calculated the overlap parameter q for small- 
world networks with various rewiring probabilities p, and 
obtained its probability distribution P(q) from MC sim- 
ulations. This distribution is shown in Fig. 0] for p = 0.1 
and 0.5 at several temperatures. At high temperature, 
P(q) shows a single peak centered at q = 0, characteris- 
tic of a paramagnetic state. The width of this peak is a 
typical finite-size effect, which should collapse to a Dirac 
5 function at q — in the thermodynamic limit N — ► oo. 
When the temperature is lowered, the distribution P(q) 
broadens, as a consequence of the appearance of an in- 
creasing number of edges displaying frustration. At still 
lower temperatures, two peaks develop in P(q), symmet- 
ric respect to p = 0, and characteristic of a spin-glass 
phase [H, [5l|, [H[ . 

Information on the "freezing" of the spins as temper- 
ature is lowered can be obtained from the evolution of 
the average value of for a given degree of disorder p. 
This average value is shown in Fig. [5] as a function of tem- 
perature for several rewiring probabilities p. It is close 
to zero at the high-temperature paramagnetic phase, and 
increases as temperature is reduced, indicatinga break of 
ergodicity associated to the spin-glass phase [H, [EH • For 
p = 0, |g| converges to unity at low temperatures, reflect- 
ing the AFM ordering present in the regular lattice. For 
increasing p, we find a decrease in the low temperature 
\q\ values, due to an increasing degree of frustration. 



B. Transition temperature 

The overlap parameter q can be used to obtain accu- 
rate values of the paramagnetic-SG transition tempera- 
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FIG. 4: Distribution of the overlap parameter q for two 
rewiring probabilities and various temperatures, as derived 
from our simulations for networks with N = 10 4 nodes, (a) 
p = 0.1 at T/J = 2.4, 2.1, and 1.7; (b) p = 0.5 at T/J = 1.7, 
1.3, and 1. 
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FIG. 5: Average of the absolute value of the overlap pa- 
rameter [(g)] for various rewiring probabilities p and several 
temperatures. Symbols are data points derived from MC sim- 
ulations on networks of size N = 10 4 . 



z 

i-i 
B 

6 



<U 
C 

s 




ture, by using the fourth-order Binder cumulant [48|, [5l| 



„(T) = if3-M% 



(6) 



This parameter can change in the interval [0,1]. One has 
gjsr = for a Gaussian distribution P(q) (high tempera- 
tures), and 5 at = 1 when \q\ — 1 (in the particular case of 
a single ground state). In general rises for decreasing 
temperature, and T c can be obtained from the crossing 
point for different network sizes N. As an example, we 
present in Fig. [5] gx(T) as a function of temperature for 
several system sizes and a rewiring probability p = 0.05. 



FIG. 6: Fourth-order Binder's cumulant gw as a function 
of temperature for small-world networks generated from 2D 
square lattices with rewiring probability p = 0.05. Symbols 
represent different system sizes N — L 2 : squares, L — 60; 
circles, L — 80; triangles, L — 100; diamonds, L — 150. Error 
bars are on the order of the symbol size. 



From the crossing point we find T c /J = 2. 175 ±0.005 for 
this value of p. 

By using this procedure, we have calculated the tran- 
sition temperature T c for several values of p, and the 
results so obtained are shown in Fig. [7J For small p, T c 
decreases linearly from the transition temperature corre- 
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FIG. 7: Transition temperature T c as a function of the 
rewiring probability p for small-world networks. Error bars 
are on the order of the symbol size. The dashed line is a guide 
to the eye. 



sponding to the AFM model on the 2D square lattice, and 
for p > 0.4 it saturates to a value of about 1.7J. Close to 
p = 0, we find a change in T c induced by the long-range 
links: T c = T^-apJ, where T c ° is the paramagnetic- AFM 
transition temperature in the square lattice and a w 2. 

It is interesting the approximately linear decrease in 
T c for increasing p up to p « 0.3. This decrease could 
be expected from the larger number of frustrated links 
appearing as p is raised. For increasing frustration, the 
paramagnetic phase is favored, and the spin glass appears 
at lower temperature (T c is reduced). This change of 
T c as a function of rewiring probability could be also 
expected from the behavior of the heat capacity shown 
in Fig. [1] For p > 0.4, small-world networks behave 
in this respect similarly to Poissonian random networks, 
in the sense that the transition temperature is roughly 
independent of p, and is close to that found for the strong- 
disorder limit p = 1 . 



C. Absence of long-range AFM ordering 

Even though all data indicate that the AFM Ising 
model in small-world networks yields a spin-glass phase 
at low temperature, one can ask if such disordered phase 
appears for any finite value of the rewiring probability. 
One could argue that some residual long-range AFM or- 
dering could be present for finite but small p values. From 
our considerations in the preceding sections, one could 
think that, for small p, the low-temperature phase still 
keeps the long-range ordering characteristic of the 2D 
regular lattice, with some defects caused by the long- 
range connections. To analyze this question, we consider 




10 10 z 10 J 10 4 10 3 10° 
Monte Carlo steps 

FIG. 8: Relaxation of the staggered magnetization on the un- 
derlying regular lattice for small- world networks with p = 0.1, 
from simulations starting from an ordered AFM configura- 
tion. The plot shows the evolution of \M 3 \ at T = 1.7 J and 
different system sizes L x L. From top to bottom: L — 40, 
60, 80, 100, 150, and 200. 



a staggered magnetization defined for the square lattice 
as usual: 



with 



M s = M A - M B 



Ma = Y, Si 



(7) 



(8) 



and similarly for Mb- Our question then refers to the 
possibility of a finite value for M s for small-world net- 
works with p > 0. 

To check this point, we have carried out simulations 
starting from an AFM ordered configuration and followed 
the evolution of M s . We prefer this procedure to directly 
calculating the low-temperature staggered magnetization 
from simulated annealing, since in this case a long-range 
ordering can be difficult to find for large nerworks, due to 
the appearance of different spin domains. Thus, we an- 
alyzed the decay of M s at temperatures lower than the 
transition temperature T c for different system sizes. In 
particular, in Fig. 8 we show the relaxation of |M S | on 
networks with p — 0.1 at a temperature T — 1.7J, well 
below the transition temperature for this rewiring prob- 
ability (T c /J = 2.10 ± 0.01). For each system size, \M S \ 
decreases from unity to reach a plateau at a finite value, 
which is clearly a finite-size effect, as seen in the figure. 
As the system size increases, such a plateau appears af- 
ter longer simulation times, and the corresponding value 
of \M S \ decreases. These results are consistent with the 
decay of the spin correlation £ (r) at temperatures below 
T c , as presented in Fig. [3] for several values of p. 
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For p < 0.1, a relaxation of M s is expected to appear 
for larger system sizes and longer simulation runs. Ev- 
erything indicates that at low T the long-range ordering 
disappears in the thermodynamic limit N — > oo for any 
p > 0. This is in line with earlier results for the FM Ising 
model on this kind of networks, in the sense that the 
paramagnetic-FM transition occurring in those systems 
changes from an Ising-type transition at p = to a mean- 
field-type one (typical of random networks) for any finite 
value of the rewiring probability p > [24J . The obser- 
vation of this mean-field character for the paramagnetic- 
FM transition requires system sizes that increase as the 
rewiring probability is lowered [30( , similarly to the decay 
of the staggered magnetization in the AFM case shown 
here. 



V. CONCLUSIONS 

The combination of disorder and frustration in the 
AFM Ising model on small-world networks gives rise to 
a spin-glass phase at low temperatures. The transition 
temperature from a high-temperature paramagnet to a 
low-temperature spin-glass phase goes down for increas- 
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